arXiv: 1112.5537v3[nucl-th] 



Nonperturbative A^A^ scattering in ^Si—^Di channels of EFT(7f) 

J i- Feng Yang^'^ 

^ Department of Physics, East China Normal University, Shanghai 200241, China 
^ Kavli Institute for Theoretical Physics China, 
Chinese Academy of Sciences, Beijing 100190, China 

The closed-form solutions of A''A'^ scattering matrices for the coupled channels '^Si—'^Di in EFT(7f) 
^SJ , with the potentials truncated at order (^(Q^yj) are presented in details with the nonperturbative 

divergences parametrized in a general manner. The stringent constraints imposed by the closed- 
■ form of the T matrices are analyzed with the help of the underlying theory perspective. The 

, realization or implementation of subtractions and the manifestation of EFT power counting rules 

(-H ■ in nonperturbative regimes are shown to differ from perturbative cases and give rise to the concept 

r-{ ■ of scenario that is intrinsic of nonperturbative renormalization of EFT. Then a number of scenarios 

of the EFT description of A'^A'^ scattering are compared with PSA data in details employing the 
effective range expansion. The problems and rationalities of various scenarios for describing realistic 
A'^A'^ scattering with large scattering lengths are presented, showing that it is favorable to proceed in 
a scenario with conventional EFT couplings and reasonable prescriptions of renormalization. The 
natural appearance of and the roles played by the fine-tuning in terms of EFT couplings are also 
demonstrated in several places and interpreted in the underlying theory perspective. In addition, 
some of the approaches appearing in recent literature are also addressed in the perspective adopted 
here. Many technical details are also given for completeness of our presentation. 
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INTRODUCTION 



<^ • In the past two decades, the NN systems and nuclear forces have been intensively studied using effective 

\ field theory method since Weinberg's proposalfll, Q and the pioneering works in Refs.jl, which provides 

ff^ ■ nuclear physics with field theoretical foundations in terms of symmetries that characterize low-energy QCD. 

■ For more about the exciting achievements and progresses in this area, we refer to the review articles, e.g., 
. However, it is theoretically fair to say that a few stumbling blocks are still in the way towards claiming 

^SJ ■ the complete establishment of the field theoretical foundations for nuclear physics. The renormalization of 

TV TV scattering is one of the intriguing issues and has attracted many authors' attention. For a comprehensive 
account of the literature on this issue, we refer to the review articles in Refs.0-Q. 

Earlier impetus to this issue was given by the discussions of the problems of Weinberg's power counting 
and the treatments of the renormalization of NN sector in Ref s.ll(]| - [l2| . Since then, a number of approaches 
have been put forward, ameliorated and discussed [T^ - [l8l [20l466l l *. accompanied with some controversies 
^ . . and debates. As summarized in Ref.j^, in short of complete agreements, there are now two main choices 

■ in renormalizing the iViV sector: (1) One insists on standard subtraction algorithm of infinities through 
expanding around some leading component of the NN potential that is first treated nonperturbatively; (2) 
The other insists on the nonperturbative treatment of the full potential up to the order of truncation, with 
the renormalization implemented following the method advocated by Lepage |67|. 

As noted by many authors El]) the problem is originated from the nonperturbative nature of the low- 
energy NN scattering: The conventional algorithm for subtraction is only established within perturbative 
contexts, not guaranteed to work beyond perturbative contexts at all. Of course, the ultimate goals and 
basic principles of renormalization, which is now best encoded in the philosophy underlying effective field 
theories, should not vary with contexts. However, its concrete implementation could be context-dependent, 
it is reasonable to anticipate that not all the patterns, contents and scenarios of renormalization could have 
been fully embodied, or manifested themselves or even foreseeable in the standard perturbative algorithms. 
Therefore, the roles played by EFT power counting rules, the implementation of subtractions and the associ- 
ated wisdoms should all be reexamined, ameliorated and adapted in order to work out a fully coherent field 
theoretical framework or algorithm in nonperturbative regimes. The present status about the NN sector 
in EFT approach implies that we are still in short of a satisfactory framework in nonperturbative contexts, 
but not the failure of the effective field theory approach. In this regard, the choices made in literature could 



* We apologize for being unable to give a comprehensive list of all the literature on this issue. 
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all be understood as various degrees of efforts or steps towards the full realization of renormalization in 
nonpcrturbative contexts or regimes, and Lepage's proposal could be interpreted as the first conceptual shift 
in this direction [67|. 

Some of the concerns mentioned above have already been touched upon somehow in our previous 
works [46i IZtI . [s^ . IsoI . [giI [63| . In this report, we wish to present a study of the NN scattering in the 
coupled channels ^Si—^Di with a more coherent account of the concerns raised above within the realm 
of EFT(;zf') with the potential truncated at order 0{p'^), the extension of the results to higher orders and 
other channels is a straightforward matter. (Recently, contact corrections to NN potentials are constructed 
beyond static approximation beginning from next-to-next-to-leading order, they are important for studying 
chiral extrapolation, but not consequential for phenomenology [6^.) It is well known that the renormalization 
of EFT(7f) has been settled in conventional means with impressive achievements (see, e.g., [1, H, [Slj). Here 
it is revisited due to the availability of rigorous solutions (i.e., closed- form T matrices) to the Lippmann- 
Schwinger (LS) equations and the feasibility of a general parametrization of renormalization prescription, 
which will provide useful insights into the issue of renormalization in nonperturbative regimes. That is, the 
closed-form T matrices serve as theoretical laboratories for exploring the crucial notions and contents of 
renormalization that are intrinsic of nonperturbative regimes and not foreseeable in perturbative regimes. 
We must stress that our analysis below does not defy perturbative like approaches at all, but highlights that 
some nonperturbative components must be incorporated one way or the other. We should also note that the 
formulation and results presented below are also applicable in many non-relativistic systems such as cold 
atoms and molecules that are dominated by short range interactions. 

As far as phenomenology is concerned, the contact potential solutions will reproduce the effective range 
theory. However, the functional dependence of the closed-form T matrices upon EFT couplings and renor- 
malization prescriptions are indeed informative of nonperturbative renormalization, a source of insights that 
is elusive in the formulation of effective range theory. Another point that is noteworthy is that an EFT upper 
scale is actually physical within the effective field theories framework as it corresponds to the real threshold 
in EFT expansion, which should NOT be vacuously taken as a regulator scale or cutoff to be removed in the 
standard renormalization algorithm of perturbation theory where no finite upper scale could be reasonably 
posited. For example, in EFT(7f), the upper scale should be of the order of pion mass Ajyj ~ m-^. For 
pionfull EFT without heavy mesons, the upper scale is naturally set by the scale heavy mesons, say, the 
mass of p meson, Aj-^-j ^ nip. The same arguments apply to other low-energy effective field theories that are 
supposed to work below specific thresholds. 

This report is organized as follows: Section [H] is devoted to the detailed formulation of the problem and 
the rigorous solutions to Lippmann-Schwinger equations (LSE) for ^Si—^Di. The renormalization of the 
integrals involved is discussed and parametrized in nonperturbative regimes. In sect. IIIII we explore the 
properties of the closed-form T matrices that show why perturbative renormalization algorithm and the 
associated notions cease to apply directly and what must be done in a nonperturbative realization or imple- 
mentation of renormalization. Also presented are the scenario notion of EFT that arises in nonperturbative 
regimes and the nonperturbative running couplings in simpler cases. Section IIVI will be devoted to a com- 
prehensive analysis of the phenomenological behaviors of the closed-form T matrices embodied in various 
EFT scenarios. The problems and rationalities for some important scenarios are analyzed with the help of 
effective range expansion and fine-tuning, with the constraints implied by fine-tunings being explored and 
interpreted. Some discussions and a summary of our report will be given in sect. |Vl 



Since the NN system scatter in coupled channels only for the triplet states, let us start with a standard 
parametrization for the on-shell partial wave S matrix for the triplet states with total angular momentum j: 



II. RIGOROUS SOLUTIONS IN EFT(7f) 



A. Preliminary Setups 




[cos 2ej (p)] exp[2i(5j^i (p)] , 
i[sm2ej{p)] exp[i((5ji;^(p) -f 



Sj±ij±i{p) 
Sj±i,j^.i{p) 



(1) 
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where and Cjip) denote the phase shifts and the mixing angle, p denotes the on-shell momentum for 

a nucleon. Under this parametrization, the on-shell T matrix introduced as S = I — ^^^T with M being the 
nucleon mass. According to Weinberg, the nucleon-nucleon potential could be systematically constructed or 
calculated using xPT,lj through counting the powers in terms oi p^ / h? or rn^/h? with A being the upper 
scale for the EFT under consideration (A ~ 0.5 GeV). Then the ofF-shell T matrices for the coupled channels 
could be found through the solving LSE's, 

Tiq\q;E)^Y{q\q)+ Jviq',k)Goik;E)T{k,q;E), Go{k; E) ^ -k^M ' 

where E^ = E + ie with E being the center of mass energy and the bold-faced letters mean that they are 
2x2 matrices in the partial wave representation: 

V ^ ( ) = v^. (3) 

Note that matrix multiplication is understood in the convolution terms of the above LSE's, namely, in 
VGoT. 

As mentioned in the introduction, we wish to work with EFT(;7f) where potentials become contact ones 
that allow us to transform the integral equations into algebraic ones following Ref. .14j to investigate some 
intriguing issues from nonperturbative perspectives, unlike in many literature where KSW expansion [lol. ITTj 
is used, see, e.g., Ref.fsij and the review articles in Refs.jl, Below we will present all the necessary 
notations and details to make our presentations self-contained. First we introduce the following notations 
for the ■^5'i— ■'Z?! channels: 

-y ( \ rp ( Tsd \ I A\ 



^ Vrfs Vdd J ' \ Tds Tdd ^ 

with the meanings of the subscripts being self-evident. Then the four individual LSE's read: 

Tss = Vss + VssGaTss + VsdGoTds, (5) 

Tsd = Vsd + VssGoTsd + VsdGaTddi (6) 

Tds = Vds + VdsGoTss + VddGoTds, (7) 

Tdd = Vdd + VdsGoTsd + VddGoTdd- (8) 
The formal solutions to these coupled LSE's are easy to find: 

Tss^{l-VssGo)-'Vss, (9) 

Tsd = (1 - VssGo)-'Vsdil - GoVdd)-' = (1 - VssGo)-'Vsd{l - GoVdd)-\ (10) 

Tds = (1 - VddGo)~^Vdsil - GoVss)-' = (1 - VddGo)-'Vdsil - GoVss)-\ (H) 

Tdd = (1 - VddGo)-^Vdd, (12) 

with 

Vss = Vss + VsdGoil - VddGo)-'Vds, (13) 

Vdd = Vdd + VdsGoil - VssGoy'Vsd. (14) 

These formal expressions, which contain UV divergences, must be so renormalized through proper means 
that the nonperturbative or closed- form solutions be preserved. Using the trick introduced in Ref.[l3|, we 

could obtain closed-form and hence rigorous solutions in case of contact potentials, where the divergent or 
renormalized parameters could be readily isolated and generally parametrized [47| . 



B. Algebraic LSE's 

To this end, we introduce the following matrices for Vxy and Txy{x, y — s,d) 



Vxy{q,q')^U^(q^)KyU{q'^), (15) 
Txy{q,q';E) = U^{q^)Txy{E)U{q'^), (16) 
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C'oiss 


C2:ss 


C'4;ss 


Ass = ^ 


C'2:ss 


C4:ss 







C'4:ss 








The lowest order 


cases 


(A = 


0,2) 



with [q^) = (1, g^, g"*, • • ■ ) being a row vector [n dimensional) in terms of external momentum q, A and r 
are the corresponding n x n matrices furnished by the contact couplings C...^ at a given order of potential 
truncation in EFT framework. Note that some of the elements of the matrices [Xxy] vanish at any given 
order due to truncation. This is an important structure that proves to be consequential, see Sect. IIII Al 

For example, at truncation order C'(Q^y-)) or A = 4, we have U'^{q^) — (1, ^ q"^), and the following 3x3 
matrices, 

0\ /O 

Asd = I C2:sd Ci-sd — A , Xdd = Ci^-dd | . (17) 

Ci,sd 0/ \ 

ipler. Evidently, the contact couplings [C...] 
Cn;...ICQ ^ Aj~," in the naive EFT power counting scheme, with A(y') being upper scale for EFT(7f). 

Using such notations, the Eas. (j5l6l7l8p can be reduced to the following algebraic equations p^. l47j. 

Tss = Ass + Xss'I{E)Tss + Xsd^{E)Tds, (18) 
Tsd = Xsd + Xsd1{E)Tdd + Xss'I(E)Tsd, (19) 
Tds = Xds + XdsT{E)Tss + XddT{E)Tds , (20) 
Tdd = Xdd + XddI{E)Tdd + Xds1{E)Tsd, (21) 

with 

where the energy dependence of the r's are self-evident and henceforth omitted. The matrix I(_E) is furnished 
with the integrals arising from the convolution with Gq, prescription-dependent at this stage. With such 
kind of algebraic parametrization, all the ill defined integrals can be isolated and parametrized as elements 
of the matrix T{E) which could be defined in any prescription, see below. This yields great convenience for 
discussing the important issues associated with power counting schemes and renormalization prescriptions. 
Moreover, the solutions to the algebraic equations are easy to obtain, which are structurally the same as the 
T matrices given in Eas. (|5l6l7l8p : 

Tss - (l-Ass2:(i;))"'Ass, (23) 

Tsd = (1 - AssX(^))"'Asd(l - I{E)Xddr^ = (1 - XssI{E))-^Xsd{l - I{E)Xddr\ (24) 

Tds = {l~XddI{E))-^Xds{l-T{E)Xss)-^ = {\-XddT{E))-^Xds{l-AE)Xss)-\ (25) 

Tdd^[l-~XddX{E)y^Xdd, (26) 

with 

Ass = Ass + XsdI{E){l - XddI{E))-^Xds] (27) 

Xdd = Xdd + Xds'I{E){l - XssAE))-^Xsd- (28) 

To obtain the T matrices, one only needs to perform the matrix multiplication U'^tU . 

As a matter of fact, we could put the above results into more succinct form if we introduce the following 
'super' matrix notations for Txy, Xxy'- 

then the coupled LSE's could be put into the following one line short equation, 

T = X + )a{E)T, (30) 



t In this report, it is tacitly assumed that these couplings are energy-independent as energy dependence in potentials can be 
removed using unitary transformations [23l] . 
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whose solution reads, 

T={l->X{E)y^\. (31) 

We have verified that Eg. ([31]) do reproduce the solutions given in Eqs. (j23l24|251 [26|) using the formulae given 
in Appendix A. 



C. Parametrization of the matrix I{E) 

The renormalization of the T matrices for NN scattering within EFT(;f') now boils down to the renor- 
malization of the matrix T{E) to be realized or implemented within nonperturbative context. To proceed, 
we will need a general parametrization of this matrix, of which a generic element formally reads 

" (32) 



(2n)^ E-k^/M 
Then, such an integral can be parametrized as follows: 



I„ = -W + 51 •^2/+ip'^"~'\ lo^ Jo +i^P, p=VmE, (33) 



Aw 
1=1 

with J... being prescription-dependent parameters (usually constants) at this stage. For example, in the hard 
cutoff (A) scheme, we have 

M f p A+p\ M A2'+i 

''°=2^(,^-2^"A^j' = -2;^27TI- (''^ 



In dimensional regularization, such integral reads, 

Jo = 0, J2/+1 - 0. (35) 



we have 



M 

Jo = -^f^, J21+1 = 0. (36) 
47r 



Actually, there is a simple and natural method for computing and parametrizing the integral in Eq. p2p 
basing on an elaborate implementation of the effective field theory strategy |69l472| : First we differentiate 
it with respect to the energy or the on-shell momentum squared (parameters external to the loop integral) 
for sufficient times to arrive at a convergent integral and carry out the integration. Then we integrate back 
indefinitely, arriving at a definite expression in terms of energy or on-shell momentum plus an ambiguous 
polynomial as a general parametrization of prescriptions to be specified or fixed later. Specifically, this is to 
solve the following differential equations 

(9^2) ^n = -» P , (37) 

the solution to this equation takes exactly the form given in Ea. (l33p with Jq and [t/2m+i5™ > 0] being the 
arbitrary constants arising from the indefinite integration with respect to jp' . Of course, the best way to 
understand the values of J... is to view them as decoupling effects of the underlying structures [4?! I69l472| | . or 
as a general parametrization of such effects. 

With the parametrization of Eq. ([5^ , the matrix X(_E) could generically be written in the following succinct 
form 

A 

X{E) = -X^U{j?)U^(j?) -f 'hi+i^Ui (38) 

1=1 
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after introducing the following symmetric matrices 



1 r ,,mm^it)] 

p'Jo dt 
Jo dt 



At/i = ^/ dt^^^^^^^^^^, (39) 
AC/,+i = 4/ dt^^^^^, l>2. (40) 



The concrete expressions for AUi at order A = 4 are listed in Appendix B. 

D. Closed-form T matrices 

As mentioned above, the closed-form T matrices could be readily obtained by sandwiching the r's in 
Ea. (|23l24l25l26l) or (|5T|) between the row and column vectors [/-^(g^) and U{q''^). Putting q and q' on-shell, 
one could readily get the on-shell matrices. After some algebra, we could find following the closed-form 
on-shell T matrices at order A = 4: 

' .r. + ^±M;4, (41, 



1 ^ Mo + lo {Vo + Mip^) + IjVip'' ^ 1 
Tsdip) T^sdP^ Tds{p) 



(42) 



1 AAq+ZoPq 

^0 + 777-nF^FrT3i' (43) 



Tdd{p) {Ni+IoV^)p 
N{D„^Vl^ + V^Na. (44) 

We note that all the parameters [A/^.., I?...] are real polynomials in terms of couplings [C...], [J2m+i,n^ > 0] 
and on-shell momentum p, which are all independent of the complex parameter Tq, for detailed expressions, 
see Appendix C. 

As a matter of fact, the functional forms of the T matrices in terms of [N...,'D..] and lo given above hold 
at any truncation order. To see this, let us start with the examination of the unitarity relation in terms of 
the matrix T furnished with the above closed-form on-shell T matrices. Inverting the matrix T, we find 

^ - \ -V^d/i-D^p'), To + -DoKV^p^) ) ■ ^^^> 
With this inverse form, it is ready to see that the on-shell unitarity is fulfilled, 

T-i_(T-i)t = T-i-(Tt)-i = ^I, (46) 

27r 

with I denoting the 2x2 unit matrix. As a matter of fact, this on-shell unitarity can be rigorously shown 
to hold at any given order of truncation for the coupled channels, for details, C.f. Appendix D. Using this 
unitarity relation as a starting point, we could also establish that the inverse form of the on-shell super 
matrix T must take the following form at any given order of potential truncation^ 

rp-l ^-T -I , ( J^ssfDss, -■N'sd/'Dsd \ /.-^ 

''^\-Msdlf>sd,Mdd/'Ddd)' ^ ' 

and [A(.., !>...] must be Io"iii*iependent polynomials in terms of [C...], [J„, n > 0] and p^, again, see Appendix 
D for a detailed proof of this corollary. 

Interestingly, Eg. dTTl) also implies that the functional dependence of the Tg^^T^^ upon Ig and [N...,'D..] 
as exhibited in Eqs. (l41l42l43l) holds at any truncation order. To see this point, it suffices to invert the 
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in (|T7|) . then one finds 

J-=To + ^°+^4, (48) 

¥7, ^ ^ ' ^ ' 

J-sd Vsd 

— =Iq + — (50) 

No = AfssNddi>^sd - ^fsd'DssVdd, Nl = NssVddVla^ 

■Do = NddVesVla. i>i = VddVssVla. T^sd = MsdVsdVssVdd (51) 
which are equivalent to Eas. (|41l42l43p . Then, (H^ follows immediately as a byproduct: 

?k±3^ . 1. (52) 
NiVo 

We note in passing that a small part of the foregoing contents have been given in cursory presentations 
in our previous works JIJ, [60l IgiI [63| . Equipped with the closed- form T matrices with their prescription 
dependence being parametrized in a general manner, we could better explore the issues around the renor- 
malization of the nonperturbative T matrices. Of course, one could also explore more physical issues with 
such closed-form T matrices. We will consider some of them in the following sections. 

III. CLOSED-FORM AND RENORMALIZATION 

In this section, we explore the closed-form T matrices obtained above to gain useful insights into the issue 
of renormalization of EFT(^) in nonperturbative regime. 



A. Isolation of lo and intrinsic mismatch 

In Ea. (j45p or (j47p . there is a very crucial structure worth emphasis: The complex parameter Iq is always 
'isolated' from the couplings [C] and the parameters [J2m+i,'7i > 0] in the inverse T^-*-. This is also true 
in an uncoupled channel for L-wave 

T,-^ = ^=Io + ^p-^^ (53) 

with Nj, , Dj, being polynomials in terms of [C...], [J2m+ii rn > 0] and , for a proof in the ^5*0 channel we 
refer to [43 • The cases of the uncoupled D channels can be verified using the above T matrices through 
letting Vsd = 0, so V^d = 0, T^d - T^, - 0, and 

This crucial structure has been addressed in our previous worksdJilleSllll. To make our presentation 
here self-contained, we further elaborate on the structure as below: (1) The 'isolation' of Xq is universally true 
in any channel and at any order, making Tq unique in EFT(7f); (2) The functional shape of the closed-form T 
matrices with respect to the on-shell momentum p is physical, hence prescription independent. Consequently, 
the 'isolation' position makes Jq become RG invariant in the closed-form T matrices, except the trivial cases 
of the 5*- wave T matrices at the lowest order Tg^ = Jq + i^j^ + ^o-s ^^^^ leads to the KSW running of 
Co;s[I3|; (3) The uniqueness and RG invariance of Jq also imply that Jq could only depend on physical 
scales, namely, the upper scale A(y) and the typical scale Qi^t/) in EFT(;f): 

M 
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Thus, Jo is no longer an ordinary renormalization scale, but a fundamental quantity to be determined 
through physical inputs in EFT(7f). 

Conventionally, UV divergences appear in the local part of the vertices with tree level couplings, and local 
counterterms can be constructed from such vertices or couplings to remove such divergences. In EFT(^), 
this means the divergences in the matrix I{E) must pair up or 'match' with the contact couplings so that 
counterterms could be constructed from these couplings. Unfortunately, such matching fails in the closed- 
form T matrices, the first example is the parameter Iq or Jq noted above. 

Actually, beyond the first two orders of potential truncation (i.e., A > 4), more and more parameters in 
I{E) become 'unmatched' in the closed-form T matrices. For example, at order A = 4, this could be simply 
seen from the following relations between some coefficients in the expansion in the factors Vq, Aq and 
Afi (see Appendix C): 

Sq-aJs = -Si;3 = -Jii^o;3 = viah- (55) 

Since any ratios constructed using such coefficients, say Si-z/S^-i or fo;3/j^i;2, should be physical or RG 
invariant, so J3 = — (5i;3/(5o;4 = —vo-.s/^ia is also constrained to be physical or RG invariant. In other 
words, J3 is also 'unmatched' with EFT couplings in the closed-form T matrices. One could also verify this 
using our previous result about the ^S'o channel where Dig^.'j,/Nis^.^2 = ~J3^ with [Dig^.r^, Nig^.2] being the 
coefficients of the highest p^-power terms in [Dis^^Nig^] respectively [47|. 

A more transparent way to look at this issue is to invert the reduced algebraic LSE in an uncoupled 
channel[6l|. It suffices to demonstrate it with the ^Sq channel at a given order of truncation: 

T^l-^^l-nE). (56) 

As no element of I{E) is zero while all the elements in the upper left triangle block of X^g^ vanish due to 
truncation, there is an intrinsic mismatch between I{E) and ArJ^^ , i.e., a mismatch between the ill-definedness 
in I{E) and the available 'pools' for counterterms from AtJ^^ within nonperturbative regime. Moreover, a 
further mismatch exists even between the nonzero entries of I{E) and X~g^: the p dependence differ. Letting 

Ar^ij develop p dependence to match the p^ monomials in I{E) would lead to nonlocal time dependence in 
Ai or the local potential Vi g^Q and in turn ruin the EFT power counting. Putting more couplings into A 
to make all elements nonzero would simply break the EFT truncation rules and could not help to remove all 
the mismatches [61] . Therefore, conventional counterterms could not succeed in nonperturbative regime and 
subtractions must be implemented otherwise, and the direct link or connection between EFT power counting 
(for potential truncation) and nonperturbative renormalization is lost. Later on, we will analyze this issue 
from the underlying theory perspective that will yield further insights. 

B. Nonperturbative 'finiteness' 

Obviously, the closed-form T matrices impose stringent constraints that preclude renormalization from 
being realized in the standard perturbative pattern. Some of the prescription parameters must be treated as 
physical parameters (or RG invariants) rather than ordinary running parameters and should be determined 
through appropriate boundary conditions or inputs of data. Then there is a concern about the feasibility of 
imposing such boundary conditions as there seems to be infinitely many divergences in LSE. 

In this connection, we note that at a given order of EFT truncation, the number of the prescription- 
dependent parameters [J...] appearing in the closed-form T matrices is finite, or the rank of I{E) is fimte[54j]. 
Amusingly, this 'finiteness' is in fact guaranteed by the truncation whose presence becomes a virtue at this 
point. Thus, we are led to a novel notion of 'finiteness' in nonperturbative renormalization: There are finitely 
many nonperturbative divergences to be dealt with. All the divergences involved in the iteration series can 
be simply factorized into a finite number of 'irreducible' divergences that furnish a finite dimensional matrix 
1{E). It is this 'finiteness' that makes the nonperturbative renormalization of T matrices feasible in principle 
and hence extends the notion of renormalizability somehow in nonperturbative regimes. 

Technically speaking, the number of 'irreducible' divergences is controlled by the highest power of momenta 
or by the scaling dimension of the contact potential, not by the times of iteration in LSE. As the scaling 
dimension of a potential is still finite in pionfull EFT due to truncation, we feel tempted to speculate that 
similar 'finiteness' may also be true in pionfull EFT. 
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After making the above observations, it remains to subtract the 'irreducible' divergences and to fix the 
residual constants using appropriate boundary conditions, which will be addressed in next subsection. Here 
we briefiy remark on a statement that renormalization is done after cut-off independence is achieved, which 
is true in perturbation theory as any residual term can be absorbed into couplings through redefinition. In 
nonperturbative cases, however, it is difficult or even impossible to absorb all the extra finite constants into 
the couplings while keeping the closed-form T matrices intact. Thus, cut-off independence does not exactly 
mean that the resulting quantities are physically faithful in nonperturbative regimes. That is why people 
often turn to perturbative approaches in practices. 



C. Underlying theory perspective 

We now turn to the underlying theory (UT) perspective. In such a perspective, an EFT just provides a 
simple parametrization and organization of all the amplitudes that are calculable in a well-defined manner 
in the underlying theory. 

Suppose we start the calculation of a low-energy (LE) process from the a complete theory that underlies 
an EFT in consideration, then the results must be well defined. If we start the calculation with the EFT 
formulation, it means that the LE projection operation (T'le) has been performed implicitly before loop 
integrations as the EFT propagators and vertices are must first be obtained from LE projection. The 
problem is that such LE projection usually does not commute with loop integrations, thus, divergences arise 
from the wrong order of operation. At one-loop level, we define the commutator of the two operation as 
below: 

C.T. = [Ae, J dfi{l)]. (57) 

Rearranging the operations in Ea. (j57p . we have, 

J dM(OAE[/(/,---)]+C.T. [/(;,•••)], (58) 

with f{l,---) denoting the integrand for an EFT process. Obviously, the commutator in JJTp just pro- 
vides the counter-term or subtraction operation needed in EFT for removing the divergences in the integral 
/ diJ.{l)VLE [f{l, ■ • ■)] in order to recover the finite integral 'Ple [/ dfi{l)f{l, • • • )] in UT. So, counter-terms 
or subtractions are intrinsically originated from the operations of LE projection and loop integrations that 
are indispensable in any sensible construction of EFT. Accordingly, subtractions must be implemented at 
loop level in effect, a natural observation from the underlying theory perspective. In perturbative contexts, 
subtractions/counterterms could be readily realized through local operators in EFT Lagrangian. In nonper- 
turbative contexts, however, such a realization is not generally guaranteed, we are forced to adhere to the 
original loop- level subtractions/counterterms or equivalent operations [l7j. 

Thus, both the closed-form T matrices and the underlying theory perspective require that subtractions be 
performed at the level of (loop) integrals in effect through whatever means that makes sense-l-. In EFT(7f), 
such subtractions are straightforward to perform in the parametrization I(i?), resulting in residual ambiguous 
parameters (also denoted as [J...]) to be fixed through appropriate boundary conditions. This is now a well 
accepted algorithm in doing the nonperturbative renormalization of EFT, see a recent review (73|. In a sense, 
the renormalization of Schrodinger equations (or Lippmann-Schwinger equations) a la Lepage adopted in 
Refs. [20l - [25l [34 - 36] can be seen as an instance of implementing subtractions at loop level and then treating 
the parameters that could not be absorbed into coupling constants as independent parameters to be fixed 
separately. 



d/i(0/(;,' 



t In our previous work [4^, the counterterms in the closed-form T matrices have been termed as 'endogenous' counterterms. 
Therefore, 'endogenous' counterterms arc just the ones that could effectively perform subtractions at loop level in nonper- 
turbative regimes. 
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D. The notion of scenario of EFT 



Now we recapitulate the foregoing analysis on the closed-form T matrices as below: (1) Nonperturbative 
divergences intrinsically mismatch with the EFT couplings or interactions in the truncated 'space'; (2) 
Fortunately, the 'unmatched' divergences are finitely many at a given order of potential truncation, so 
the whole issue is still tractable and hence renormalizable in nonperturbative sense; (3) The 'unmatched' 
divergences have to be subtracted at loop level with the residual parameters being physical and independent 
ones [ J.^.'^'"'"''' ] ; (4) Then, a novel concept intrinsic of nonperturbative formulation, the 'scenario of EFT' 
(iSeft), naturally emerges in the parametrization of the closed- form T matrices [s^ l63j : 

5kpt = [c...(/i)]e{[j.(.'"^=^]®[j...(/i)]} 
= {[c...{p)]®[J...{^,)]}®[J^^^''^] 

= [d.':''='='] ® [j.^.-^'^^^i, (59) 

with ji being the running scale in EFT(7/). Namely, a scenario consists of EFT couplings and complementary 
parameters that only arise from loop integrals; (5) Direct link between EFT power counting and nonpertur- 
bative subtractions/counterterms is lost due to the mismatch mentioned above. EFT power counting could 
only be manifested through the renormalized EFT couplings, not directly applicable to the bare objects or 
counterterms in nonperturbative regimes [isj. 

Therefore, we anticipate that the scenario structures just delineated above must be effectively incorporated 
somehow in the EFT descriptions of the NN sector. In this regard, at least the approaches adopted in 
Refs. pol - [25l . [s^ - iSGl lss I all essentially incorporate such notion of scenario. For example, the separation scale 
A (^ 750MeV) in Ref. [5g] essentially plays the roles of the additional parameters that are complementary to 
EFT couplings in contrast to the EFT running scale /i that is of order to^. In other words, the discrimination 
between /i and A is necessary and natural in the light of EFT scenario. Of course, such notion of EFT scenario 
could not be fully appreciated in the Wilsonian RGE analysis of EFT couplings as it is focused on EFT 
couplings, not on the full scenario 'space' [62j. The indispensability of the complementary scenario parameters 
could only be seen in the 'full' (closed-for m) amp litudes or T matrices. In practice, such parameters could 
appear in various disguises, see, e.g., Refs.[20l425l l33 - [39l [Sll - fssl [ssj. or even embodied somehow in different 
specification of 'dynamical' degrees (ssllSGj. 

In our analysis performed so far no deformation or extra construction has been introduced into the EFT 
framework. All the notions and observations addressed above are merely natural consequences of the closed- 
form T matrices that could not be foreseeable in perturbative context. Hence, the original goal of EFT 
approach-providing field theoretical foundations to nuclear physics-is preserved. Throughout the whole 
analysis, the stringent constraints imposed by the closed-form T matrices are not circumvented, instead, 
they are directly confronted, exploited and finally turned into virtues in understanding the whole issue 
within field theoretical framework. 



E. Nonperturbative running couplings 

As was shown above, the closed-form T matrices require all the parameters generated from loop-level sub- 
tractions to be discriminated into two subsets: one collects the physical or RG invariant parameters that are 
complementary to couplings, the other consists of the running ones that are intertwined with the couplings. 
Hence, the couplings develop nonperturbative running behaviors, which could be easily demonstrated with 
the coupled channels ^Si—^Di at order A = 2, where Jq is the only physical complementary parameter: 

Ml = 0, AAo = (1 - Co;ss Js)', 'Dl = - C2;.d(l - C2;ssh), Vo = So;0 + <5o;ip', with 

Then, according to Sect. IIIlXl the following combinations arc RG invariants besides Jo: 

<^0;0 ^0:1 a "^sd 
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then the T matrices can be physically parametrized as below: 



1 _ 1 

~ -^0 H ; 5 — T- o2 4' 



J_ _ ^ 1 + (gp + a2 P^)To 
J_ ^ l + (ao + a2p')Xo-Xg/?V 



= ^0 + ^^55;;4 , (63) 

(64) 



From these physical parameters, we could find the running couplings that absorb the running parameters J3 
and J5: 

Co,. = Jf I'^U ' ^^'^ 

C2,ss = (^1 - , C2;.d = |, (66) 



^ = ^l + a^h-p^Jl (67) 
These nonperturbative running couplings possess both IR and UV fixed points in literal sense 



= 0, = 0, C^-} = 0. (69) 

Like the ^Sq case[47j, these nonperturbative running couplings all vanish in UV regions as they are only 
low-energy couplings that will naturally be dominated and superseded by the couplings that govern the 
interactions in UV regions. Actually, the EFT description breaks down before the running scale goes to the 
UV end as these nonperturbative running couplings blow up already at a finite J3: 



Q2 ± ^Jaj + 4/3^ 

Jz;± ■ (70) 

This property also corroborates the fact that this EFT only makes sense below a physical upper scale. There 
is no point in letting the running scale or the subtraction point take very large values. In pionfull EFT, of 
course, it is an extremely challenging task to correctly calculate the contributions from all the intermediate 
states, especially the sophisticated suppression of all the modes beyond the upper scale. 

IV. VARIOUS SCENARIOS OF EFTC/) 

A. High- and low-energy behaviors 

First, let us entertain ourselves with some interesting estimates about the high- and low-energy on-shell 
behaviors of the closed-form T matrices obtained above, which is straightforward to see from their transparent 
p dependence. 

From Appendix C, one could easily read off the following high-energy or UV on-shell behaviors (i.e., p — > 
large) at truncation order A = 4: 

^^j^ + ^^p + o{p-% (71) 
Tss 47r 

^ = '^o+i^P + o(p-2), (72) 
Tsd = Tds = o{p-''). (73) 
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In terms of the parametrization defined in Eq.(IT|), this is: 

^"f (p) =(ns + ^^TT + o{p-^), ns e Z, (74) 



53^,^ ip) = (j^d + + o{p-'), Ud e Z, (75) 

e(«-)(p) =n,7r + o(p-^), £ ^. (76) 

Of course, such high-energy behaviors and their significance could not be taken seriously for the realistic 
NN scattering in terms of EFT(y) which is only a valid description well below 1 GeV. However, as a general 
theoretical formulation for non-relativistic systems, such behaviors might be physically reasonable in certain 
systems that will be studied elsewhere. The point we wish to make is that, such high-energy or UV behaviors 
are compatible with unitarity. 

Meanwhile, in the low-energy or infrared limit (p 0), we have 

^ = Jo + f^ + o{pl (77) 

Jss 00;0 



o^o-.o 



Tsd = Tds = o(p2), (79) 



and 



^35,(0) =ns7r, nsGZ, (80) 
<^3Di(0) = nnTT, no e Z, (81) 
e(0) = UETi, UE e Z. (82) 

For the realistic NN scattering, we know that ns — l,nj:i — ue — 0. Obviously, these behaviors are also 
compatible with unitarity. These low-energy or infrared limit behaviors are directly related to the effective 
range expansion (ERE) within the realm of EFT(/-). We will examine the ERE parameters in details in 
Sect. UvBl 

We note that the foregoing behaviors are obtained from the expressions of the parameters [Af...,!)...] as 
monomials of on-shell momentum squared (p^) at order A = 4, where the highest power 'w' of each monomial 
can be read off from the expressions listed in Appendix C: w-Dq = 4, ujui = tOAfo = 3, wa^^ — i^v^d — 2- Then 
it is rational to expect that the following ranking should hold at higher order of truncations: 

min{ajx)i , i^Afo } > ^v.^ ■ (83) 

Then one may be convinced through the above rankings that the above limiting behaviors should qualitatively 
persist at higher orders of truncation. 



B. Scale hierarchy and scenarios 



In the following, we examine the behaviors of the closed form T matrices obtained above in the light of EFT 
scenario and extend the preliminary analysis given in Ref.j6^. Obviously, different values of the scenario 
parameters will lead to different ratios constructed with the coefficients in the expansion of [A/"... , I?...] , hence 
different physical behaviors in terms of the closed-form T matrices. Reversely, the correspondence between 
physical behaviors and these ratios also constrains the choices for EFT scenarios. To study the physical 
behaviors of the T matrices and their phenomenological consequences, we first perform some magnitude 
order analysis of various parameters involved. Here and below, the subscript '(t^)' ^(i) ^'^'^ 
omitted to avoid heavy symbolism. 

Generically, the scenario parameters [C] and [J...] depend on the ratio e = Q/A or /i/A that stipulates 
EFT expansion due to scale hierarchy: jj, ^ Q, Q A. For the realistic NN scattering in the realm of 
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EFT(;f) where A ~ m^, a-^{^Si) ~ 36.4MeV, a'^{^So) ~ -8.3MeV, etc., we could envisage the following 
scale hierarchy: 



1 1 
e — - 



4 ■ a(3^i) 



6A, 



1 



-e^A,.... 



To proceed, we introduce the following dimensionless parameters for the couplings and J.. 



47r C2n;...(e) 
M 2"A2»+i 



47r 



J2fe+i(e), 



(84) 

(85) 
(86) 



where C2n;- - may be 0(1) multiplied by powers of '2' due to our convention. For the complementary param- 
eters is replaced with 'Q^fe+i' in 



In this report, we will consider the following three typical scenarios for simplicity: 

M 

A: £2n,.. - 0(1); J2fc+i - 0(i); Jo - ^0; 

B: C2„;... ~ i2fc+i ^ C(i); Jo ^ ^Q; 

C: C2„,.. 0(1); Wi ^ 0(1); Jo - ^A. 

47r 



(87) 
(88) 
(89) 



Obviously, scenario A is given in terms of conventional EFT power counting for couplings and a usual 
prescription for [J...], it will lead to natural ERE parameters, hence a natural scenario. Scenario B looks more 
complicated as unconventional power counting of couplings is incorporated: The couplings are unnaturally 
large. This scenario will indeed lead to unnatural scattering behaviors. In scenario C, everything is simply 
the same as in scenario A except Jo(= Re[Zo]). This is because as a physical or RG invariant parameter (see 
Ea. ((M)) ). Jo could simply be a function of the physical upper scale A only[43,|ll,|6l|. Or, due to its uniquely 
'isolated' position in the closed-form T matrices, we promote it to be a fundamental parameter that scales 
as C(^^s in scenario C! 

Actually, scenario C is 'natural' in the sense that all the scales involved are 'naturally' sized, but it 
could also lead to unnatural scattering lengths for ^-waves upon reasonable fine-tunings without using 
unconventional EFT power counting, see below. In a sense, our formulation (in the simple scenario C) 
provides a 'physical' and natural foundation to the EFT treatments (see reviewsj^-Q) that employ various 
forms of fine-tuning. Recently, it is emphasized that fine-tuning is at least partly responsible to the complexity 
in nuclear physics ^7^]. In our view, the field theoretical origin of such complexity or fine-tuning lies in the 
nonperturbative nature of EFT renormalization, or, specifically, in the nonperturbative scenario of EFT. 
Further remarks on fine-tuning will be given in Sect. IIV El 



C. T matrices in various scenarios 



In the closed-form T matrices, fine-tunings may 'naturally' occur among the scenario parameters and lead 
to some intriguing phenomena (like unnatural scattering lengths). For our purposes below, it is convenient 
to define the fine-tuning in a scenario in terms of the dimensionless couplings [c...(e)] only: 

=±l + o(e'^), a€(0,K), (90) 



|C2„;...(U 

where k denotes the power exponent of the next-to-leading order contributions to the coefficients [v..., 5..] 
(C.f. Appendix C) in e expansion. Evidently, it makes no sense to let the fine-tuning exponent a take any real 
number larger than k. This is the guiding principle for determining the choice of fine-tuning in the following 
discussions. Larger k means larger capacity for fine-tuning, less sensitivity to higher order contributions, 
and finally more credits for EFT approach in that scenario. The exponent k in a scenario could be obtained 
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through the examination of the e dependence of the T matrices. For purpose of demonstration, we wiU 
simply choose j2k+i = 1 or 



in the following calculations. 



1. Scenario A 

With the above preparations, one could find that the on-shell T matrices for the ^Si~^Di channels would 
depend on the expansion parameter e as follows after replacing /i and Q with eA, 

47r/T,, ^ ^ ^ ip ^ l + o(6=^) + ^0(l + o(e-^))+--- 

(91) 

Closer study shows that for Cq-^ss in this scenario: Ka = 3, which means a large capacity that favors stable 
fine-tuning against higher order corrections. However, this favorable capacity for fine-tuning is useless, which 
could be seen through the following simple computations in ^Si channel: 

AO(l), re - 2— — . (92) 

a co-ss A A 

That is, the naturalness of these two effective range expansion (ERE) parameters could not be altered by 
any fine-tuning in this scenario or system. Therefore, to obtain unnatural scattering length within the realm 
of EFT(7f), one must alter either the EFT power counting for couplings (like scenario B) or the scaling laws 
of the parameters [J...] (like scenario C), or both. 



2. Scenario B 



We note in passing that this scenario mimics the so- called KSW[l3 

scheme in the part of couplings within 
the realm of EFT(7f) while taking the same prescription of J... as in scenario A. Like above, in this scenario, 
we could obtain the following e dependence of the T matrices. 



4^/T,, , ip , l + o(e) + f^O(l + o(e)) 



A 5o;.. + o(e) + ^0(l + o(e)) + 



(93) 



Note that, in our convention, C2n; -- ^ scenario B. Here, one could find that Kb = 1 for cq-^ss, 

which means much less capacity for fine-tuning, an unfavorable point for this scenario. Nevertheless, an 
unnaturally sized scattering length could indeed be achieved for ^Si channel in this scenario as a result of 
the unconventional rating of the couplings: 

a \ cq-ss J A A 

However, this scenario would lead to other unnaturally large ERE parameters in S'-channels also due to 
the unconventional rating of couplings, see next subsection. As a matter of fact, it is possible to achieve 
unnatural scattering lengths without modifying the conventional power counting: To tap the capacity stored 
in the complementary parameters. To this end, we consider below the simplest case of such scenarios. 
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3. Scenario C 

In this scenario, we could find the following e dependence in the T matrices: 



(95) 

With such T matrices we find that Kc = 3 for co;ss, a large capacity for fine-tuning. More pleasantly, this 
large capacity for fine-tuning is pivotal here for producing large scattering length. To see this point, we 
perform similar computations as above and find that 

1^-a(i + ^)^A.o(.^), ...2%^^. (96) 
a \ co;ss J A A 

Here the fine-tuning co;ss = ^ 1 ^ o(e'^) with a E (0,3) is employed. For the realistic scattering, it 
suffices to choose cr = 1 so that a^^ ~ eA. Thus, in the realm of EFT(y), a natural effective range and 
an unnatural scattering length could be 'naturally' achieved in '^S'l channel in scenario C. Naturally, it is 
more instructive to compute and compare more ERE parameters across various scenarios. This will be done 
below. 



D. Effective range expansion in various scenarios and PSA data 



To proceed, we start with the ERE in L-wave defined as below: 

p^^+^ cot 6 l{p) = - - + IreP^ + 
a 2 

n=2 

In ^Si^^Di, the phase shifts are given as below in terms of the parametrization given in Eq.(IT]): 

1 _ 4Mp^ 



^I^tUp) 



l-^\TsdipW 

'l-^\Ts,{p)V 
After some algebra, we could find the following simple relations: 



p cot5s{p) 
p^ cot 5d{p) 



in 
An 



Re 



1 



I. 



dd\ 



o[p' 



(97) 

(98) 
(99) 

(100) 
(101) 



according to which, we could compute the ERE parameters up to Vi in '^Si simply using Tss, while in 
it is legitimate to compute the ERE parameters up to V2 with Tdd- The resulting expressions for these ERE 
parameters must be rational functions in terms of [z^..., J...] and Jq. For example, the scattering lengths for 
^Si and ^Di can be expressed exactly in terms of vo-fi, ^o;o, i^i;0) and Jo as below. 



1 



a(35i) 
1 



a(3Di) 

which could be further expanded in terms of e. 



An vq.q + Jo(5o;o 
'M So;o 
An i^Q.Q + JpSo-o 
M + Jo(5i;o' 



(102) 
(103) 
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The detailed e dependence of the ERE parameters is further determined by the fine-tuning of relevant 
couplings in a scenario. The results about the first two orders in e expansion are summarized in tables I and 
II, where the following two primary fine-tuning patterns (tuning I and II) for the leading coupling cq-ss are 
demonstrated as we are at this moment mainly concerned with the scattering length for ^Si: 

Timing I: cq-ss ~' — 1 — o(e) (scenario A,C), eco-ss ~ — 1 — o(e) (scenario B); 
Tuning II: co-ss — 1 — o(e^) (scenario A,C). 

Note that it does not makes sense to consider the tuning II in scenario B as it is simply not allowed there 
due to Kb — 1. Fine-tuning of higher couplings will be considered later for higher ERE parameters. We 
should also note that in order to yield a scattering length of order (eA)~^ in scenario B, one should use 
ecQ-ss ^ +1 + o(e) instead of tuning I, with the rest being essentially not affected. Thus, different scenarios 
would embody different tuning mechanisms to yield the same scattering length. Of course, they will not 
always agree with each other about the higher ERE parameters. Then for one particular system, only one 
scenario will survive experimental tests. Interestingly, as 

^i!^ = ^-M (104) 

these two scattering lengths could be both naturally sized (scenario A) or both unnaturally sized (scenario 
B and C). 



TABLE I: NaturaIness(N)/unnaturaIness(U) of ERE parameters in ^Si—^Di:ciy^ss ~ — 1 — o(e) 



ERE 


Scenario A 


Scenario B 


Scenario C 


S: A-a 

A - re 

A^ -V2 
A^ - V3 


C>(l + o(e)) 
2c2-as + o(e) 
0{l + o{e)) 
0(l + o(e)) 
0{l + o{e)) 


e~^0{l + o{e)) 
2e^C2-ss + o(e) 
e-'0{l + o{e)) 
e-^0(l + o{e)) 
e'^0{l + o{e)) 


6-^0(l + o(e)) 

2C2:ss + o(e) 

0(l + o(e)) 
0(l + o(e)) 

O(l + 0(£)) 


D-. A^ -a 

A^'' - re 

A-l • V2 


0(l + o(e)) 
20(1 + 0(6)) 
0(l + o(e)) 


e~'0{l + o{e)) 

20(1 + 0(e)) 
e-'0{l + o{e)) 


6-^0(l + o(e)) 
26-^0(1 + o(e)) 

6-^0(1 + 0(6)) 



TABLE II: NaturaIness(N)/unnaturaIness(U) of ERE parameters in ^ S\~^ Dr.co-ss ~ -1 - o(e^) 



ERE 


Scenario A 


Scenario C 


^Si: A-a 
A - re 

A"" -V2 
A^ ■ V3 

A'' -Vi 


0(l+0(6)) 

2c2;ss + o{e^) 
0(l + o(e)) 
0(l + o(e)) 
0(l + o(e)) 


.-^0(l + o(.)) 

2&2:ss + 0{e^) 

0(l + o(e^)) 
0{l + o{e')) 
0{l + o{e')) 


^Di: A" -a 
A-3 • re 

A-l ■ V2 


0(l + o(e)) 
20(l + o(e)) 
0(l + o(e)) 


e-^0{l + o{e)) 
2e-2O(l + 0(6)) 
.-*0(l + o(.)) 



For comparison, we also present here the analysis for the ^5*0 channel in table Hill where to yield a much 
larger e~^) scattering length, the tuning eco ^ — 1 + o(e) is used in scenarios B while in scenario C we 
employ the tuning cq ^ — 1 + o(e^). 
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TABLE III: Naturalness(N)/unnaturalness(U) of ERE parameters in ^So 



ERE 


Scenario A 


Scenario B 


Scenario C 


A 


a 


0(1 - 


ho(e)) 


£-^0(1 - 




£-^0(1 


+ o{e)) 


A- 


Te 


2C2 


^oie) 


26^52 4 


o{e) 


2£2 + 


o(e^) 


A3 


V2 


0(1 - 




£-iO(l^ 




0(1 + 




A^ 


«3 


0(1 - 


Ko(e)) 


£-^0(1 - 




0(1 + 




A^ 


«4 


0(1 - 


ho(e)) 


£-=^0(1 




0(1 + 





In spite that the foregoing analysis appear to be quite crude, the basic features of the three typical 
scenarios are clearly revealed: Scenario A characterizes systems with natural scattering behaviors, while the 
rest two account for unnatural systems. For TV TV scattering, to see the rationality in choosing appropriate 
prescriptions rather than modifying the canonical EFT power counting, we perform some further magnitude 
analysis of the ERE parameters in '^S'l and channels of the PSA data used in Ref . [TSj] . The results are 
given in table HVl where we denoted the upper EFT scale and scaling parameter as A and e respectively, the 
actual values used in our analysis A and e are charged pion mass m^± and j. From this table, one can see 
that the PSA data lead to small V2 in comparison to all the three schemes above. In particular, the PSA 
data give an extremely small V2 in ^5*1 channel. 



TABLE IV: ERE parameters in S-waves: PSA data. 



ERE 


'•^Si: data 


scaling 


^So: data 


scaling 


A • a 


(0.26)-^ 


£-10(l) 


-(0.06)-^ 


£-^0(1) 


A • re 


(0.81)-i 


0(1) 


(0.53)-^ 


20(1) 




(4.13)-^ 


£^0(1) 


-(1.81)-=* 


e4o(l) 




(1.53)-^ 


eio(l) 


(1.07)"^ 


0(1) 




-(1.16)-^ 


eio(l) 


-(0.92)-^ 


0(1) 



Comparing tables HI and HIl with the PSA data, it appears that the predictions for W2,W3 and U4 of S*- waves 
in scenario B are 'wrong': they are unnaturally 'large'. So there are huge 'gaps' between PSA data and 
scenario B: 

: !M ^ ,-4^ !M ^ ,-5, ^ (105) 

V2,P V3-P V4-P 

iQ ^2;B _9 V3.B _2 «4;S _3 ..„„^ 

So- , e (106) 

V2:P V^.P V4.P 

with subscript '... ;b' for scenario B while '... ;p' for PSA. In the meantime, such 'gaps' are relatively smaller 
in scenario C: 

3^,: !^^,-3, !!Mi^,-f, !M^,-f, (107) 

V2;P V3-P V4.P 

^So: ^^e". (108) 

V2;P V3-P V4-P 

Thus, scenario B seems to be disfavored in comparison with scenario C as the latter is relatively closer to 
the PSA data. In ^Sq, the agreement between scenario C and PSA data is almost complete. 

As the numbers in tables HI [II] and IIIII had been derived with simple fine-tunings of the lowest coupling 
Cq;... only, one may wonder what would happen after taking the fine-tunings of higher couplings into account. 
Looking at the actual expressions, one could find that higher ERE form factors [vk,k > 2] involve more 
higher couplings at the 'leading' order of e expansion. This implies that cancelations may occur amongst 
the couplings involved after considering appropriate fine-tuning so that the resulting ERE form factors could 
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be reduced in magnitude. As will be shown in next subsection, this could indeed happen in scenario C. 
In scenario B, however, such fine-tunings of higher couplings are hard to manipulate for two reasons: (1) 
On the one hand, there are much larger 'gaps' to be reduced, which in turn requires much larger capacity 
for fine-tuning; (2) On the other hand, there is actually little capacity for fine-tuning to be exploited as 
Kb = 1. So, this conflicting tension between larger 'gaps' and little capacity for fine-tuning means that the 
magnitude-reduction of higher ERE form factors in S-waves is practically impossible to achieve in scenario 
B. The little capacity for fine-tuning in fact stems from the unconventional rating of couplings in scenario 
B, i.e., they are larger than in scenario A or C. In this regard, scenario B, or a scenario with unconventional 
power counting of couplings, is strongly disfavored in the EFT description of NN scattering. 



TABLE V: Low energy theorems in KSW 



ERE 




scaling 


'So 


scaling 


W3(fm^) 


-0.95 
+4.6 

-25 


0(1) 


-3.3 
+17.8 
-108 


^ A3 

e-f £(1) 

^ AT 



More than a decade ago, treating pion exchanges perturbatively in NN scattering using KSW power 
counting was shown in Ref.fisl IT6| to lead to some bad low energy theorems (LET), see table [V] These 
let's are qualitatively in agreement with our rough estimates in scenario B in that the ERE form factors 
derived there are also much larger (at least one order of e~^) than the PSA data: 

5^1 : — — e ^ , — — ^ e ^» , — — ^ e ^« , 

V2:P Vs-p V4:P 

'So ■■ — e e 9 , — e 3. 

V2-P Vs-p V4.p 

Therefore, the scenarios with unconventional power counting of couplings where couplings are promoted 'up' 
in magnitude seem to be problematic choices for NN scattering in lower partial waves. 

E. Fine-tuning and hidden structures 

In this subsection, we discuss the feasibility of deriving smaller higher ERE form factors in scenario C. It 
suffices to demonstrate it with V2 in '^Si which is the smallest ERE parameter in S"- waves according to the 
PSA data as listed in table HVl 



1. A small V2{^ S\) in scenario C 
Let us first look at the detailed expression of V2 in scenario C in leading orders of e expansion: 



A^-V2^ '--^ + o{e'), (109) 



from which it is clear that for the higher couplings C2-ss, C2-sd, C4-ss and c^-ss, we still have a large capacity 
K = 3 that could also be seen from the coefficients of 'p^' and 'p''^' in the fraction part of T^^^. That means, 
the fine-tuning exponent a for these couplings also lies in the interval (0,3), just like that for cq-ss- Then, 
a small V2 of order 'o(e'^)' would be resulted if the couplings conspire to essentially cancel out against each 
other on the right hand side of Ea. (|109p up to order o(e'^). Considering V2 alone here, this could be achieved 
in a number of ways. For example, the following choices in combination with tuning I for cq-ss could result 
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in a V2 of the required size: 



C2-SS ^ l~ o(e^) 
C2-sd o(e2) 
Ci-ss ~ 1 + o(e) - o(e2) 

li-ss 2 + 2o(e) - 2o(e2) 

co-ss ^ -I- o(e) 
A^W2 ^ o(e^). 



(110) 



(111) 



Of course, other ERE parameters may impose further fine-tuning requirements for the couphngs, which 
would effectively constrain the freedom in the choices of fine-tuning for all the couplings involved. 



2. Combtned constraints 



The foregoing analysis shows that, to conform with the empirical sizes or values of the ERE form factors, 
the EFT couplings must be further constrained through fine-tuning. In other words, the couplings must be 
correlated somehow in the context of fine-tuning. Below we would like to provide another instance of deriving 
constraints or correlations for the couplings: The combining considerations of the fine-tuning requirements 
from the scattering lengths in ^5*0 and ^5*1 in scenario C. From our results given in the foregoing subsections, 
it suffices to consider the couplings at the lowest order. 

To proceed, we adopt the following decomposition of contact NN potential at the leading order of 
truncation [zi, : 

VnIi = Cc + {Tl-T2)Cl,c + {cTl-CT2)Ct 

+ (ti •T2)(cri •,T2)C/;t. (112) 

Projecting it into partial waves, we have the following lowest order couplings for ^Sq and ^Si, 

CiSo = Cc + Cl-e-iiCt + Cl-t), (113) 

as, = Cc + Ct-3{Ci.e + Ci.t). (114) 
Then from tables I and III, the ERE analysis requires following equalities in scenario C 

^Ci5o = -l + o{e'), (115) 

^Cs, = -l-o{e). (116) 

Now, it is evident that, the four lowest couplings [Cc, Ct, Ci;c, Cj-t] are correlated with each other in the 
leading order of e expansion: 

Ct{0) = Ci-AO), (117) 
CciO) - 2C/;e(0) + 3C/;t(0)-l. (118) 

In the EFT approach to NN scattering, these relations follow as corollaries in scenario C. They may be tested 
using other methods like lattice computation. Furthermore, if the assumption that ^j^Cc{0) — -^^Cf (0) = 1 
is reasonable, then the above relations imply that the fourth coupling C/ t is suppressed by at least one order 
of e: Ci-tiO) = 0, Ci-t = i^o(e'^), a > 1. 



3. Interpretation 

Now we make some attempts at the understanding and perhaps the interpretation of the above deductions 
concerning fine-tuning. We have shown that fine-tuning in combination with phenomenological requirements 
can lead us to arrive at quite some 'orders' or 'structures' hidden in the EFT couplings. These hidden 
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TABLE VI: ERE in variants of scenario B: eco-.ss ~ — 1 — o(e) 



Variant 


j2fe+l = 1 


j2fc+l = 


:• jfc+2 

J2fc+1 — £ 


S:A-a 




£-^0(1 


£-^0(1) 


A ■ re 




2e2c2;,,C'(l) 


2£2C2;..0(1) 


A3. ^2 




6-^0(1) 


£-10(l) 


A'^ . i;3 


£-^0(1) 


£-^0(1) 


£-^0(1) 


A^.t;4 


£-30(1) 


£-30(1) 


£-30(l) 


D : A^ a 


e-*C'(l) 


£-''0(l) 


£-*0(l) 


A-3 . re 


20(1) 


20(1) 


20(1) 


A-^V2 


£-30(1) 


£-30(l) 


£-30(l) 



'regularities', should they be correct or trustworthy, must come from the 'order' or 'symmetry' contents of 
the underlying theory. 

Let us elaborate. It is natural to expect that the contact couplings of EFT(/-) should be proportional 
to Ijm"", thus such constraints must reflect the structures of the pion exchanges expanded in terms of 
momenta against pion mass, hence the structures of broken chiral symmetry intrinsic of chiral perturbation 
theory or QCD as the underlying theory. As a matter of fact, this situation could be generalized as below: 
Each coupling in an EFT is a simplified projected version of the amplitude computable in the underlying 
theory. These amplitudes must be constrained with fundamental 'symmetries', thus 'correlated' with each 
other. Translated into the EFT language, then, we end up with the EFT scenario parameters fit in certain 
fine-tuning patterns that are reminiscent of the correlations of the underlying theory amplitudes. Thus, EFT 
and the associated fine-tunings (sound ones, of course) are nothing else but a coarse or simplified realization 
or reflection of the delicate and yet regular 'structures' in the underlying theory. Not knowing the details of 
the underlying theory, we have to constrain the scenario parameters with empirical data or other physical 
inputs, thus the fine-tuning in the EFT approach to nuclear forces is in fact an indispensable step in the 
course of imposing boundary conditions. In particular, in the scenario C considered here, fine-tuning is a 
'fine' 'component' of EFT calculations, naturally driven and soundly supported by empirical data. 

V. DISCUSSIONS AND SUMMARY 

In the foregoing analysis, we have worked with the uniform prescription for [J2fc+i]. One might wonder 
what would happen if alternative prescriptions are adopted. As these parameters should not be characterized 
by scales larger than Q or /z, then possible choices could only be characterized by smaller [J2fe+i], such as 
j2fc+i ~ f-''^^ or even j^k+i — 0. With such variants, one could find that the above conclusions for scenarios 
A, B and C essentially remain intact. For scenarios B and C, we summarize the ERE form factors calculated 
in the variant prescriptions for [J2fc+i] in tables |Vl] and IVlT] respectively, where one could find that the three 
versions of scenario B and C essentially yield the same ERE form factors. Closer studies would tell us that 
all properties presented above remain valid across the three versions of scenario B and C. In particular, the 
leading order expression of V2i^Si) in the three versions of scenario C is identically the same, hence the same 
fine-tuning works in the three versions of scenario C that could reduce the size of this ERE form factor from 
0(1) down to o(e3). 

The results given in table I VI ll reveal that the PDS like prescription [j... = 0) essentially produces the same 
effects for ERE as the usual one {j... — 1) does, provided Jq is separately determined as a physical scale. That 
means, in the PDS like prescription of [j...], conventional EFT couplings are compatible with large S-waves 
scattering lengths in the closed-form T matrices provided Jo is properly sized. This reminds us again that 
Jo is a fundamental parameter in the EFT description of NN scattering. As this parameter arises from 
the rescattering diagram, thus it must be a crucial parameter associated with rescattering. At this point, 
we refrain ourselves from attempting at further interpretation of this parameter in terms of nuclear physics 
before more phcnomcnologically oriented investigations are performed, especially when pion exchanges are 
present. 

Below, we wish to make some attempts at understanding the various approaches proposed and/or adopted 
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TABLE VII: ERE in variants of scenario C: co;ss ~ — 1 — o(e) 



Variant 


J2fc+1 = 1 


j'ik+i = 


l = 6^"+^ 


S -.A-a 




e-^0{l) 


6-^0(1) 


A • re 


2C2;ssO{l) 


2C2;ssOil) 


2C2;ssOil) 


A3.«2 


Oil) 


Oil) 


Oil) 


A'' • «3 


Oil) 


Oil) 


Oil) 


A^ ■«4 


Oil) 


Oil) 


Oil) 


D -.A^ - a 


e-'0{l) 


e-^Oil) 


e-^Oil) 


A-3 • re 


2e^^Oil) 




2e-^Oil) 


A-i ■ 


e-'^0{l) 


£-^0(1) 


6-20(1) 



in literature in the light of EFT scenario discussed so far. Roughly speaking, a vast number of papers 
employed a finite cut-off of various sorts to remove the divergences. According to our foregoing discussions, 
the rationality of such treatments as field theoretical framework could be assessed in the light of EFT 
scenario: The finite cut-offs essentially play the roles of the complementary parameters, whose presences 
implement the loop-level subtraction or introduce the 'endogenous' counterterms effectively. In other words, 
as noted in Sect. IIIIEl the EFT description usually breaks down at a finite value of the running scale, thus 
there is essentially no point to include the intermediate modes at scales much higher than the upper scale 
of EFT [HO, HI]- Otherwise, such treatments could not be seen as genuine field theoretical ones due to the 
modeling 'elements' built-in with a finite cut-off that carries no physical essence in EFT perspective. 

On other hand, there is also a vast number of perturbation like treatments in literature that could also 
be understood in terms of EFT scenario (see, e.g., Refs. [iol |49|. [stI [ssj ) as the complementary parameters 
are also introduced or incorporated there in various disguises. Typically, they are either introduced as 
separation scales and/or 'allocated' somehow in the additionally introduced couplings that are promoted 
according to certain modification of EFT power counting. As such parameters could not be simply absorbed 
into conventional couplings or taken as ordinary running parameters, the scales involved in these parameters 
are practically fitted to some values near the upper scale of EFT or appropriate empirical data. 

In a sense, the two main choices mentioned above seem to converge somehow to the EFT scenario structure 
explicated in this report: the conventional EFT description in terms of power counting only is augmented by 
complementary parameters in various forms, e.g., finite separation (or cut-off) scales fitted somehow and/or 
additional couplings that are promoted according to modified EFT power counting. Roughly speaking, 
keeping the regularization-related scales finite is just an alternative way of incorporating the complementary 
parameters into EFT description, hence an indirect realization of nonperturbative scenario of EFT renor- 
malization. Of course, one thing that is common in all these literature is that the 'leading' components must 
be treated nonperturbatively. Moreover, fine-tuning is also at work to various extent in all these literature. 

Similarly, the approaches adopted in Refs.|30l l5ll - l53| could also be interpreted as incorporating the notions 
such as loop-level subtraction and complementary parameters that could not be readily identified as or 
absorbed into conventional couplings. Of course, the scenario structures beyond conventional couplings 
might also be realized by incorporating unconventional degrees and the associated couplings [5g|. In 
the contexts beyond few-body systems, it was noted in Ref.f?^ that the perturbative like treatments of 
pion-exchanges still has some problems to fix, while the other main choice seems more efficient [7l-[9l [79| . 

In summary, the closed-form T matrices for NN scattering in the coupled channels "^Si—^Di were presented 
in a general parametrization of divergent integrals within the realm of EFT( /f ) . The tight constraints imposed 
by these closed-form T matrices led us to the nonperturbative scenario of EFT renormalization: There are 
intrinsic mismatches between the EFT couplings and divergent integrals that can be factorized into finitely 
many nonperturbative divergences to be treated, then subtractions must be performed at loop-level through 
whatever means which give rise to some complementary parameters that are constrained to be physical or 
RG invariant. In addition, the couplings develop nonperturbative running behaviors. Then several typical 
scenarios were examined and analyzed using effective range expansion. It was shown that the scenarios with 
unconventional couplings suffer from some shortcomings and also seem to be disfavored by PSA data, while a 
simple scenario with conventional couplings is indeed compatible with realistic NN scattering behaviors after 
appreciating the pivotal role played by the RG invariant parameter Jq and the fine-tuning of couplings. A 
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heuristic understanding of most literature on NN scattering was also presented in terms of the fundamental 
contents of the nonperturbative scenario delineated here. We stress again that, our investigation has been 
performed in a general way that is applicable to any consistent EFT dominated by contact or short-distance 
interactions, not just valid in chiral effective theories of low-energy NN scattering. More phenomenological 
works such as descriptions of phase shifts and mixing angles using our closed-form T matrices for pionless 
EFT will be done in the near future [8fl]. The relations between the closed- form T matrices and 'perturbative' 
ones in KSW expansion will also be examined in future works [8l|. 
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Appendix A 

Suppose Xab{A^ B — 1,2) are four n x n matrices, then the super matrix 

Xll Xi2 \ 



X = 

\ X21 X22 



has a formal inverse of the following form. 





X 


-1 


[(2L 


-;)n> (^-;)i2] 

"^)2i ' ^)22 / 






11 


= Xii- 


-'^12-'f22^"'^2l) 






12 


= {X21- 


- X22^12^"'^ll) ^' 






21 


= (^12- 


- XiiX2j^X22) ^, 






22 


= (-^^22 " 


- X2iX^j^Xi2) ^. 



(119) 

This is valid when each of the four blocks is non-singular matrix. When Xab 7^ B) are singular, we have, 

{X_ = (-'^12-'^22^^21 - -'^11) ^Xi2X2^, 

)21 

or 



^ ^)91 — -'^22^^2l(-'^12-'^22^^21 - -'^11) \ (120) 



{X_ ^^Y2 ~ X^^Xi2{X2lX^^Xi2 — X22) 

{X'^)21 = (^22-^21^n^l2)"'^21^n . (121) 



Furthermore, if three of the sub matrices are singular, then X ^ does not exist at all. 
For 1 - }X[E) in Sect. IH we have 



-1 _ / A^ss, ^ss^sdK-dd 

K-dd^ds^ss, K^dd 



{i-\i{E)r = 

with 

^xx — (1 ^xx-^) 1 ^xx — (1 ^xx-^) : 

where Xxx is defined in Eq. ()27|28p and x — s,d. 
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Appendix B 

At order of A = 4, we have 
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Af/i = 


1 




P^ 


, AU2 = 
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'1 


, AC/3 = 










, A[/4 = 














P" 


P'J 






p' 
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Appendix C 

Introducing the following parametrization of [Af...,!)...]: 

3 2 4 3 2 

AAo = ^ ,yO;jP^' ,J^1 = Y1 Vl-jp"' ,Vo = Y, ^0;jP^' ' ^1 = E ^^-^ ''^sd = Yl ^-d'^jP^' ' (122) 
j=0 j=0 j=0 j=0 j=0 

the concrete expressions of the coefficients [i'..., S...] read as follows, 

Z^0;0 = (1 ^ C2:ss^3 ^ C'4;ss Js)^ — (C'4;ss + C^-ddji^ — Ca^-^ssJ^Y Jb ~ C'o;ss C'4;ss j| + 2(1 — C4:ss>/5) 
X [(C2;.ssC4;drf ^ C2;sdC4;sd) 'h — (Ci-ssCiiSS + C4:sdC'4:sd) "Z?] ^3 + [(2C2;.ss C2;.sdC'4:sd — C2-gg 
xC'4:dd ^ C'2:sd'54;ss)^5 + 2C4:sd (C'2;ss C'4;,,d — C2:s(iC'4;,,s) Jy — C^-ssiC^.gs + C%,.j}Jq\J1 
+ {Ci-ssCi.dd — C'4;sd){(l ~ Ci-^ssJbf' JI + Co-ssJ^Jb + 2C4;ss J3 J5 J7(l " C'4;ss>/5) 

+Cl;sd'^l{J^J^ - J7) + (^1,15^3 J5J9}; (123) 
I'O;! = { — C'4;ss(l — C'4;ss'^5)^ — C'4;dd(l — C2-ssJif' — "^{Oi-gs + C2-sdCi-sdJi){^ — C2;ssJ3) + 2(C|.gg 
+2C4;ssC'4;dd — C4,-sdCi-sd)Jh — C2;sd^4:;ssJi + 2[(C2;ssC'4;sd — C2;sdC'4;ss)C4;sd + 2(C2;sd 
xC'4;sd — C2-ssC4;dd)Ci-ss]J3J5 — ^Ci-ddCl-sgJl + "^Ci-ssCi-sdCi-sdiJt + -'S'^?) + C'4;ss(C4;ss 
~^4;sd)'^3'^7 + (C'4;ssC'4;dd — C'4;sd)[^0;ss'^3 + 2(1 — Ci-ssJhf Jh + ^Ci-ssJsJ? + C'4;ss'^3('^3 

X J9 - 3 J5 J7) + C|.,rf( J3 Jg - J5 J7) J3]} J3; (124) 

l'0;2 = {C'4;ss + (C4;ss(^4;dd — C'4;sdC'4;sd)(2 — 2C2;ss'^3 " 3C4;ss'/5) + (C4;ssC4;sd — C4;s(iC4;ss) 

X (2C2:sd^3 + Ci-^sd'h) + (C'4:ssC'4;dd — Cj.sdJK^ — Ci-gsJhf — {Cl-ss + ^4;sd)'^3J7]} Js'i (125) 
J^0;3 = {2Ci-ssCi-sdC4;sd — C^.^gC^^dd — CA-ssCXgd)^^'^ (126) 
^1;0 = C'4;(id(l — C2;ssJ3 ~ C'4;ss>-^5)^ + 2C2:s(iC'4;sd(l — C'4;ss^5)^i ^ (2C2:s.sC2;.sdC'4;.sd — C'lisd 
xC'4;ss) j| + (Cl-^ji — C4-ddC4;ss)[Co;ssJ3 + (1 " Ci-ssJbf'Jb + 2C4;ss(l — Ci-ssJblJsJl 

+{Clss + Clsd)JlM-, (127) 

= {2(C4;sciC4:sd ^ C'4:s.sC4;(M) (1 — C2-,ss'h ^ C'4;.s,s>/5) + '2'C2-^sd{Ci-^sdCA;ss — C4;ssC'4;sd)<^3 

+ (C'4:sd - C'4;ddC'4;ss)[(l - .-/5C'4;ss)^ + (C|;ss + C'4;sd)'^3'/7]}'/3; (128) 
i^l;2 = (C'|.^sC4;(id + C'4;ssC'4;sd ^ 2C4:ssC4:sdC'4;s(i)>/3 ; (129) 
^O-fi = [Co-ss + {Cl-ss + C'l;^^) J9][l - (Qiss + C4^-dd)Jb\ + (Cfiss + C^2;sd)'^5 + '2{C2-ssCi-ss 

+C2-sd04-sd)J7 + {C2-sdC4-sd — C2-ssCA-dd){C2;ssJt + 2C4;ss J5 J7) + {C2-ssCi-,sd — (^2;sd 
xC'4;ss)(C2;sdJ5 + '^Ci-sdJ^Jl) + {Cl.ssCi-ss + "^Ci-ssCi^sdCi-sd + C%sdCi-dd)J7 + {C^-ss 

^Ci-dd - Cl.sd)[Co;ssJ'i + {Ci-ss + (^4,sd){JiJ9 - J5J7)]; (130) 
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So;l — 2C2;ss + (C^sd ^ ^2;ss)'^3 + (C'4;ss — C4-dd)[Co;ssJ3 + {Cl-gs + ^l-scdJ^J^ — Cl-g^J^Jj] 

+'^C2;sdCi:sdJ5 + '^{C2:sdC i;sd — C2-ssC4-dd)[J5 + C4;ss'/3>^7] + (C'4:ss + C'4:s(i + 2C4-ssC4:ss 
+2C4-^sdC4-sd)J7 + 2{C2issC4-sd — C2:sdC'4;ss )C'4;sd ("/s " J3J7) + [2C4;ssC'4:sdC'4:sd — Cl.,,g 
x{C4-ss + C4-dd)]J5J7 + {Ci-sd ~ C'4;ssC4;dd)[2C4;ss J5 + (C|.^s — C'f.^^;) ( J3 J7 — ^5)]^?; (131) 

^0;2 = 2C4;ss + (74:55 + (C4:ss + C4:dd)[(l ^ C'4;ss>/5)^ ^ 1 ^ 2C2:ss^3 + {C^i-gd ^ ^4:ss)'^3>^7] 

+2(C4;sd + C4-^sd)C2;sdJz + [C'4;s(J ^ C'4;ss + 2C4;sciC'4;sd)]«/5 + C2;sd (C'2;s(iC'4;ss ~ C'2;ssC'4;sd) 
X j| + (C2;ssC'4;d(i — C2-,sdC4-^sd){C2;ssJl + 2C4;ss>/3>^5) ^ ^l^sd^i-ssJ'^ + (W-^^ — C'4;ssC'4;d(i) 
X{(1 - C4;.. J5)V5 + 2(1 - C4;ssh)C4:sshJ7 + [Cq:ss + + Q^d) Jg] J3 } ; (132) 

^0:3 = {C\;sd ^ ^i:ss + 2 (C4;ss ^4;^^ + C4-^sdC4-^sd) + 2 (C4;ss C'4;c;d — C'4;sdC'4:sd) ((^2:8^ •/3 + C4-^ssJb) 
+2(C4;s(/C'4;ss — C'4;ssC'4;s£i)C2;sd^3 ~ {C4;ssC4-^dd ~ C'|.s^)[(l — C4:ss>/5)^ ~ (C'4;ss + ^4;s£i) 
XJ3J7]}J3; (133) 

'^0:4= (C4;ss^4;dd + C'4;ss(74;s(i ~ 2C'4;ssC'4;sdC'4;s£i)>/3 ; (134) 

(^liO — [Co;ss + (C4:ss + Cl.,.^)Jci\C4-^dd " C2:sd + [^2;ss^4;d(i + C^.^d^^^ss ~ '2C2;ssC2;sdC4-^sd]Jb 
+2(C2:ssC'4;dd — C2;sdC4;sd)C4-^ssJ7 + (C'4;ss C'4;d(i " W.^^) [(74.3^ J7 — Co;ss>/5 " ((74.3^ 

+C2,,)J5J9]; (135) 

(^l;! = 2[C2:ssC'4;(Jd — C2;sd{C4-sd + C'4;sd)] + (C'4;ss C'4;dd — Cl.^^;) [Co;ss >/3 + 2C4:ss>/7 + (^^l-g^ 
+^4;sd)'^3>^9 + (C'l-sd — C|.^g) J5 J7] + {C2;sdC4;sd ~ C2;ssC4-^dd)C2-,ssJ3 + (C2:s(iC'4:ss 
— C2;ssC'4;sd)(2C4;s(i>/5 " C2-,sdJ'i) + [((^f.^s — C\.gd)^i;dd — "^G 4-^ssC 4-^sdC 4;sd]J7] (136) 

1^1:2 = {C4-^ssG4-^dd — Gl._gd)[{^ ~ C'4;ss Js)^ — {G'4-_ss + ^4;sd)'^3>^7] + "^iC 4;ssC 4-^dd — C4-^sdC4;sd) 
x(l — C2-ssJ3) — Cl.gd + 2(C4;ssC'4;sd ^ C4-sdC4-^ss)C2-,sdJ3 + (C'4;ssC|.g^ — C^f.^^ 

xC4;dd)J5; (137) 
(51:3= — ((7|.g^C'4;dd + C'4;ssC'|.s^ — 2C'4;ssC'4;sc;C'4;sd)^3; (138) 
<5sd;0 = C'2;sd(l — C2-ssJ3 — C4-^ssJb) + [Co;ss + [Cl-ss + (^f.^,^) J9]C'4;s(i>/3 + (C'4:ssC'4;sd + C4-^sdC4-^dd) 

x(l — C4-^ssJb)J7 + {C2;sdC4-^sd — C2;ssC 4-dd)C 4-sdJ3J7 + [C2;ssC'4;sd " C2;sdC4-^ss + {C\.sd 

-~C4-ddC4-ss)C4-sdJ7][{^ ~ C4-ss'h)J5 + C4;ss J3J7]; (139) 
<5sd;l = C'4;sd + C'4:sd — [C2-^ssC4-^sd + C2;sdC'4:ss + (C4;sd + C'4;dd)C'4;ss C4;sd >/7] ^^3 " [C'4;ss C'4;sci 

+C4;ssC'4;sd(l — C'4;ss>/5) + C4;ssC4;sd]</5; (140) 
5sd:2 = —C4-ssC4-sdJ3] (141) 



and 



VQ;3 + !^1;2 J'3 = 0, 5o-4 + ^l;3>/3 = 0, ^0;4 = Vl;2- (142) 



Appendix D 



In this appendix, we present the rigorous proof of the lo-independence of the fractional parts of T^-*- in 
■^^i— '^-Di. To this end, we need to prove the unitarity of T first. 

Using the super matrix notations introduced in Sect. IIIBl the algebraic LSE's read, 

T = A + AIr. (143) 

Due to the EFT truncation, the enlarged matrix A and hence r are singular, their inverse do not exist. Thus, 
we need to prove the following unitarity relation in terms of r 



.Mp 



UU^ 

uu^ 



(144) 
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or more specifically, 

f Tss Tsd\ _ ( Tss Tsd \ ^ -^^P ( T^s T^d \ ( UU'^ \ '^ss ^sd \ 
\ Tds Tdd / V Tds Tdd / 27r I Tds Tdd \ UU'^ / I Tds Tdd j 



We start with the Eq. (|143l) . and its hermitian conjugate, 

rt = A + TtlU, (146) 

where the transpose symmetry of A has been used. Now, multiplying Eq. (|143p from the right by It} and 
multiplying Eq. (|146p from the left by tI \ we could find that. 



tIt} = AIrt + Alrlrt ^ tIt} = -t} + (1 - AZ)"Vt, 
rTrt = tIX + TlT^t}\ rlV^ = -r + r(l - 1} X)-^ . 



(147) 
(148) 



Noting that 



we finally have 



(1-AI)"V^ = (1 - AI)"^ A(l -lU)-i =t(1 -ZU)-\ (149) 



liZ! - T)l! = t} -T. (150) 
In similar fashion, we could also prove that 

z^d^ -Sx = z^ -z- (151) 

This ends the proof for the unitarity in terms of t. We have also verified this relation using the solutions of 
T... in terms of A... given in Eqs. (j23l24l25l26p . 
Now sandwiching the unitarity between 

\ C/ / \ u I 



we find 



U^TdsU U^TddU j \U^rdsU U^TddU J \Tds Tdd J \Tds Tdd J 



\ t 

sd \ 



.Mp / U^T,, U^Tsd ] UU^ W f/^T,, U^T. 

' 2n [ U^Tds U^Tdd )[ UU^ )[ U^Tds U^Tdd ) 
27r I U^TdsU U^TddU i I U^TdsU U^TddU 27r ^ ^ 



Then inverting T again, we arrive at 



T-i-(Tt)-i=i^I. (153) 
27r 



These are the unitarity relation in terms of T matrices. Since Jq always go with i^^, we conclude 

T-i=IoI + AR, (154) 
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with AR being Zo-independent. In fact, each element of AR takes the form of ^ with Af and T) being 

polynomials in terms of [C...], [J2m+i] and p^. The arguments go as follows: It is easy to see that T is a 
matrix made of rational functions in terms of [C...], [J2t?x+i] and p^, thus must also be such kind of 
matrix. As the only complex parameter Iq has been isolated, the rest must be a real rational matrix in 
terms of [C...], [J2m+i] and p^. Q.E.D. 
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